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In stars and planets, magnetic fields are believed to originate from the motion of electrically 
conducting fluids in their interior, through a process known as the dynamo mechanism. In this letter, 
optimisation of the velocity field for magnetic energy growth enables us to simultaneously address 
two fundamental questions of dynamo theory: "Which velocity field leads to the most magnetic 
energy growth?" and "How large does the velocity need to be relative to magnetic diffusion?" . Here 
the full space of continuous solenoidal velocity fields in a periodic box is considered. Measuring the 
strength of the flow with the root-mean-square amplitude, an optimal velocity field is shown to exist, 
but without limitation on the strain rate, it is locally optimal only. Measuring the flow in terms 
of its associated dissipation leads to the identification of a single optimal at the critical magnetic 
Reynolds number necessary for a dynamo. This is found to be only 15% higher than that necessary 
for transient growth of the magnetic field, implying that stationary velocity fields are sufficient at 
low magnetic Reynolds number. 

PACS numbers: 47.20.-k,95.30.Qd,47.54.-r 
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The continuous stretching and folding of magnetic field 
lines by a velocity field is considered to be the main mech- 
anism generating magnetic fields in stars, planets and 
the interstellar media [1]. This magnetic dynamo mech- 
anism must counter magnetic diffusion, which occurs on 
a time scale that can be estimated by L 2 /\, where L 
is the length scale of the system and A is the magnetic 
diffusivity. Together with a scale for the velocity U, the 
relative growth versus diffusion can be estimated with 
the magnetic Reynolds number, Rm = LU/X. With- 
out complete knowledge of the interior flows of astro- 
physical bodies, theoretical studies have considered many 
parametrised velocity fields in several geometries, includ- 
ing the Ponomorenko Q, Roberts [J, ABC 0, Q and 
Dudley and James flows [6j, over the last 40 years. Sys- 
tematic searches continue to seek the best dynamo pos- 
sible dynamos within the parameter space, e.g. [7j. 

The small length scale of the laboratory compared 
to astrophysical bodies implies rapid diffusion, and the 
mechanism is therefore difficult to reproduce. Neverthe- 
less, this is an exciting era where laboratory experiments 
have begun to realise this process The experi- 

ments vary greatly in geometry, but in order to be suc- 
cessful, all seek to optimise the flow conditions necessary 
to realise magnetic energy growth. 

In this letter, optimisation is shown to be possible 
without need for the specification of a parametrised set 
of acceptable flows. This enables a lower bound on the 
magnetic Reynolds number to be identified for a dynamo. 
Optimising for both long- and short-term growth reveals 
that the dynamo mechanism need not rely on transient 
growth effects to work well, and that a stationary velocity 
field is sufficient at low Rm. 

Throughout this letter, the length scale is taken to be 
L = L x j(2ix\ so that the scaled box has length 2ir in 



each direction, and the following notations are used: 
v 2 = v-v, (a) = -^JadV, \\a\\ 2 = (a 2 )^, (1) 



where V is the volume of the box, so that || • ||2 is equiv- 
alent to the root-mean-square value. To begin with, the 
velocity scale is taken to be U = ||^||2- 

A variational optimisation method is used to find the 
velocity field u that maximises the growth in B after a 
period of time T. This method has recently proven useful 
in the study of the growth of disturbances in shear flows 
11 , 12| . Consider the objective function 
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The first term on the right-hand side is to be maximised. 
The remaining terms are constraints, including Lagrange 
multipliers A^, 11^ = Ui(x) and T = T(x). These terms 
are enforced to be zero. After applying variational deriva- 
tives it may be written that 
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"ind." represents the induction equation, as it appears in 
([2]), and "adj." is set to zero giving the adjoint equation 

-W = (Vxr)xu + -L v 2 r . (7) 

Km 

In deriving these expressions it is necessary to lift deriva- 
tives off the variations, for example 
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where the product rule and Gauss' Theorem have been 
used. For the first of the final two terms, the integral 
over the closed surface vanishes for case of the periodic 
box. For the second, it is quite beautiful that the La- 
grange multipliers themselves provide projection func- 
tions — these will be used to ensure that u and Bo are 
solenoidal. 

For a given u and £?o, both solenoidal and normalised, 
timestepping the induction equation to give Bt ensures 
that the penultimate term on the right-hand side of ([3]) 
is zero. Then dC/SBo in (|3j) and (j6j) is set to zero 
with the compatibility condition Tt = 2Bt- Given IV, 
timestepping the adjoint backwards sets the last term 
in (|3]) to zero and provides IV Now all quantities are 
known to calculate ascent directions for C given by (j4]) 
and ([5]) . New fields that lead to an increased C are given 
by u := u+e (SC/Su) and Bo •= Bo+e (5£/5Bo), where 
e is a small scalar value. The new fields are projected 
onto the space of solenoidal functions by considering the 
divergence of the update, which defines the projections 
functions IV After this, the are then chosen such that 
the new fields have unit norm, and all constraint terms 
in (j2]) are then zero. 

The value of e is adjusted according to whether or not 
consecutive updates appear to be pointing in a similar 
direction. Note also that B needs to be known for all in- 
termediate times during the backwards integration of T. 
This could require significant computer memory. Instead 
B may be saved at 'checkpoints' and re-integrated for- 
wards when need, involving only 50% extra work overall. 
Spatial discretisation used in the time-stepping code is 
via a triple- Fourier expansion. Nonlinear terms are eval- 
uated pseudospectrally on a grid with at least 36 points 
in each direction. This resolution was found to be more 
than sufficient for the majority of calculations, where Rm 
is very low. 

Figure [T] shows the result of calculations starting from 
several random initial u and Bo at Rm = 1. All con- 
verge to the same optimal, and the error, measured by 
(((SC/Su) 2 } + ((SC/SBo) 2 })^ drops by 5 orders of mag- 
nitude during the calculation. A sure test that the end 
state is an optimal is to verify that perturbations lead 
to a decrease in the growth. The sensitivity to pertur- 
bations of the optimal is shown in the inset to Figure [TJ 
where the growth rate a of the field B is estimated at 
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FIG. 1. Several initial conditions converge to the same opti- 
mal; Rm = 1, T — 1. Inset: Sensitivity of the optimal u op 
at Rm = 1.737 to perturbations u p measured by the growth 



rate, a in units H^l^/^. Here u 
velocity fields have unit norm. 
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FIG. 2. Magnetic energy for converged optima at Rm =1.5, 
1.737 and 2.0. By time T the signal is dominated by the lead- 
ing eigenmode. Inset: Divergence for Rm = 2.2; quantities 
from bottom to top are cr, H^l^ = IMI2, ||<*>||oo and H^Hoo, 



time T = 7. At this low Rm the optimal is apparently 
quite robust. 

Increasing Rm, the optimal u op is found to change little 
until at Rm = 1.737 the growth rate a is zero. Figure [2] 
shows the dependence of (B 2 ) with time. At these low 
Rm only modest transient growth is observed. By time T 
the energy is steady and Bt is dominated by the leading 
eigenmode. In this non-dimensionalisation the diffusion 
time L 2 /A is equal to t = Rm in units L/||tx||2. Helicity 
is plotted in Figure [3] to give a sense of the geometry 
and symmetry of the flow. The magnetic field, however, 
is concentrated at a subset of the stagnation points in 
the flow; see [13[ for plots where the dynamo mechanism 
appears to be similar. 

Although the optimal can be traced a little beyond 
Rm = 2, it becomes clear that its basin of optimal- 
ity rapidly becomes vanishingly small, as not all op- 
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FIG. 3. Optimal at Rm = 1.737. (a) Isosurfaces of posi- 
tive and negative helicity u • cj (yellow, green) and stagnation 
points (blue). (6) Isosurfaces of (pink). 

timisations converge. While the scaled velocity satis- 
fies 1 1 1 1 2 = 1, the strain rate S = (2 Sij Sij)* , where 
Sij = \{pi Uj + dj Ui), is unlimited. Large S is known to 
lead to large growth. The inset to Figure [2] shows an ex- 
ample where the method deviates from the local optimal 
at Rm =2.2. Improved resolution does not remedy the 
problem, it only permits even finer structures to form. 

A more appropriate measure of the flow is therefore 
necessary, here taken to be the velocity scale U = L | \S\ {2- 
Whilst S and the vorticity uo = \uj \ are not equal locally, it 
can be shown that ||5'||2 = 1Mb- The latter is somewhat 
more convenient to work with in the variational method. 
The scaled velocity field satisfies ||c*;||2 = 1, so that opti- 
misation can be considered to be over the space of fields 
with equal viscous dissipation or input power driving the 
flow, neglecting any feedback on the flow from the mag- 
netic field. This scaling leads to the magnetic Reynolds 
number 

Rm. = ^» (9) 
In a similar context, Backus [14[ derived a lower limit 



for dynamos in a sphere in terms of the analogous mag- 
netic Reynolds number except involving ||S , || 00 rather 
than I I2. Although not used further here, higher norms 
of S and uj are observed to behave similarly; see the inset 
to Figure [2j 

The optimisation is only slightly altered, where now 
C=(B 2 T ) - A!((a; 2 )-1) - ... (10) 
and the new update is given by 

— - = / Bx(Vxr)dt-2AiVxw + Vni. (11) 

Similar to before, the scalar Ai is chosen such that for the 
new u one has (uj 2 ) = 1. The extra complexity of the 
update leads to a linear approximation that is valid over 
a shorter range, and the number of iterations necessary 
for convergence is typically order 1000, compared with 
order 100 before. 

Upon optimisation with the vorticity scaling, almost 
the same optimal was found at Rm^c = 2.48. Its corre- 
sponding Rm = 1.75, now an observed quantity, is only 
slightly higher than found previously, where is was op- 
timised for Rm. Starting from 20 initial conditions at 
Rm uc , all converged to the same optimal. Figure |4] com- 
pares growth rates for the optimal found at Rm = 1.737, 
where u is fixed and Rm^ is varied, with the optimal that 
could now be tracked up to Rm^ = 100. For this range 
of Rm^ the fixed velocity field is competitive with the 
optimised state, which changes relatively little. From 
Rm^ = 2 AS to 100 it remains an optimal with only a 
relative change of 8%. The inset shows that the op- 
timal remains fairly robust to perturbations. Beyond 
Rm u = 100, convergence was found to be possible using 
a higher norm ||cj||4, but this incurs further iterations for 
convergence. Optimisation at higher Rm becomes com- 
putationally expensive — in addition to increased spatial 
resolution and more iterations, larger target times T are 
necessary to pass longer transients. 

Time dependence of the velocity field is an important 
factor that could affect magnetic energy growth. In order 
for energy growth to be enhanced, the changing velocity 
field must exploit transient energy growth, to beat the 
mean of the energy growths associated with each velocity 
field considered separately. At low i?m, however, this 
affect has barely been observed. Setting T to a small 
value permits calculation of velocity and magnetic fields 
that lead to the largest initial magnetic energy growth. 
This is shown in Figure [5] using T = 0.05 and 20 initial 
conditions. Two optima were identified, the lowest sets 
the energy stability bound at Rm U g = 2.12. 

In summary, a minimum magnetic Reynolds number 
is found for a kinematic dynamo at Rm uc = 2.48, by 
optimisation over the space of velocity fields with equal 
||cj||2 or associated viscous dissipation. Starting from 20 
random initial conditions at this Rm^, it was the only 
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FIG. 4. Growth rates, a in units ||c*;||2, for the optimal at 
Rm = 1.737 (red, solid), where the velocity is held fixed. 
This is bounded above by the optimisation at each (blue, 
dashed); Rm^c = 2.48. Inset: Growth rates for perturbations 
to the optimal at = 50; disturbances as defined for Fig- 
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FIG. 5. Development of energy following optimised initial 
growth; Rm^g — 2.12. The local optimal for which growth is 
initially zero defines Rm^g (red, solid), but it leads to rapid 
decay immediately following the initial 'burst'. A second op- 
timal (blue, dashed) decays monotonically but more slowly. 
This second optimal is connected to that for long-term growth 
for Rmuj > Rm^c- 

optimal found. This velocity field is very close to a lo- 
cal optimal at Rm = 1.737, in the space of fields with 
\\u\\2 = 1; it is a local optimal only, as the strain rate 
is unlimited. A lower bound for instantaneous magnetic 
energy growth is found at Rm^g = 2.12. This is not far 
below Rrriujc = 2.48, the lower bound found for long- 
term growth induced by a steady velocity field. The 
small difference between Rm^g and Rm uc leaves little 
room for improvement to be gained by time-dependent 
velocity fields. Transient growth will, however, be in- 
creasingly important at large Rm. While this velocity 
field already appears to be a fast dynamo, with a growth 
rate a = 0.357 ||<*>||2 = 0.527 | \u\ \2/L, optimisation at 



larger Rm will improve this further. 

As shown here, the variational method simultaneously 
determines a lower bound on the Reynolds number for a 
dynamo, whilst identifying velocity fields that maximise 
magnetic energy growth. It is a rare occasion that we 
can put a numerical figure to an important parameter. In 
geometries closer to experiments, tricky boundary condi- 
tions on the magnetic field for spheres and cylinders pose 
technical challenges, but it will be worthwhile overcom- 
ing them. It will be interesting to see if velocity fields 
realisable in experiment are close, or can be made closer, 
to optimal. That the growth rate is observed to change 
little for perturbations about the optimal, is promising 
for the dynamo's robustness. It will also be of interest to 
further assess the mechanism that optimises growth, and 
optimisation at higher Rm may identify an alternative 
optimals and dynamo mechanisms to that observed here. 

The author thanks Eun-Jin Kim for helpful discussions 
and Chris Pringle for an introduction to the variational 
method. 
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